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ABSTRACT 


The computer code that was developed as a result of 
the analysis of Part I of this report, VPI-AERO-063, is herein 
described and, documented. Included are programs input, ex- 
pected output , description of mode of operation and other in- 


formation useful in running the program. 
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NOMENCLATURE 


metric coefficient of the € coordinate 

metric coefficient of the w coordinate 

(heh,,) “ah, /2, curvature of the — coordinate line 
(heh) "ah, /2E, curvature of the w coordinate line 
local body radius in stagnation point polar coordinates 
surface distance along & coordinate 

inviscid velocity component in € direction 

inviscid velocity component in w direction 


inviscid velocity component in &€ direction at the body 
surface 


w= W, at stagnation point and on symmetry planes, 
W= Ue elsewhere 


inviscid velocity component in w direction at the body 
surface 


distance along stagnation point axis 
W/U, 

W,/U, 

W/W 

(h,U,) E 9U,/ 9 

ee a E 3U,/3w 

(h, We) 
(heW.) 


orthogonal surface coordinate in longitudinal direction 


1 


E dW,/ dw 


1 
E OW, /d€ 


source density 


meridional angle in stagnation point coordinate system 


orthogonal surface coordinate in transverse direction 


INTRODUCTION 


A computer program has been developed for solving 
the three-dimensional boundary-layer flow on blunt bodies 


using a finite difference procedure. The program is capable 


of calculating laminar or turbulent flows and accounts for 


the effects of surface curvature. The equations solved, 
finite difference procedure used and results of test calcu- 
lations with the code are presented in Ref. l. 

The boundary-layer program is not a stand alone 
program, but requires the generation of a data tape which 
describes the body geometry to be solved and the surface in- 
viscid pressure distribution over the body. This report is 
divided into two sections, the first describes program TAPGEN 
which is capable of generating the required data tape for the 
boundary-layer code, and program ICBL3D which actually solves 
the boundary-layer equations for the body of interest. 

All of the programs have been coded in FORTRAN IV 
and written for the IBM 370/158 computer with the JES2 opera- 


ting system. 


SECTION I 


PROGRAM TAPGEN 


In this section a brief description of the computer 
programs used to generate the inviscid data is presented to 


enable the user to execute the programs easily. The input 


data preparation and the outputs available have been described 


in detail. Finally the JCL required to execute the programs 
on IBM/370 Model 158 with JES2 operating system is listed. 
The necessary computer programs are HESS, AOAT, BLOT, 
DERVAT and INVTAP. Program HESS generates the body geometry 
and inviscid velocity distribution for two onset flows: one 
in x-direction and the other in y-direction. These two flows 
are combined by program AOAT to give the inviscid velocity 
distribution at the given angle of attack. Programs HESS and 
AOAT may be replaced by program HESS3D for calculating the 
inviscid flow over bodies with non-axisymmetric shapes. Pro- 
gram BLOT reads this velocity distribution and computes the 


position of surface coordinates, metric coefficients and the 


inviscid velocities in surface coordinates. All this informa- 


tion is read by program DERVAT which computes the required 
derivatives. Program INVTAP fits appropriate Fourier series 
in »; computes the Fourier coefficients and writes them on a 
tape as required by the boundary layer programs. These pro- 
grams taken together represent program TAPGEN. 
1. Body Geometry and Inviscid Velocities 

The boundary-layer equations have been written in 
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an orthogonal curvilinear surface coordinate system with its 
origin fixed at the stagnation point. The boundary-layer 

computations need the body geometry, metric coefficients and 
the inviscid velocity distribution. These are computed sepa- 
rately and a brief description of these computations follows. 

The surface of the body is defined in the cartesian 
coordinate system shown in Fig. 1 . These coordinates form 
the input to the computer program which computes the inviscid 
velocity distribution at the given angle of attack. This com- 
puter program is capable of handling axisymmetric bodies at 
angle of attack and three-dimensional bodies at zero-angle of 
attack. 

The inviscid flow velocities obtained above have been 
transformed into a orthogonal curvilinear surface coordinate 
system in which boundary-layer computations were made. The 
position of the surface coordinate system and the metric coef- 
ficients have been determined as required by the boundary-layer 
program. The mathematical details of these computations can 
be found in Ref. 2 

Since the position (x, r, ¢) of the surface coordinates 
and the position of the inviscid velocity would not be the 
same in general, the inviscid velocities at the surface 
coordinates have been obtained by double interpolation technique 
in (x, >). Then the results are written on a scratch tape unit 
for every £=constant. The quantities on tape are $4, hey ho» > 


s, U,, Uy and x. Since the boundary-layer program needs the first 
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derivatives with respect to &, of h;, hy: U,, Uy» x, r and the 
second derivative of r, they have been computed using the 
five-point Lagrangian formula and written on the same tape. 
These data are fitted with appropriate Fourier series in 9 
and the Fourier coefficients have been written on another 
tape. It is from this tape that the boundary-layer program 
reads and obtains the necessary information about the geometry 
and inviscid velocity distribution of the body. 
Immediately following is a description of the input 
to the inviscid program package along with sample input/output 
and a description of the Job Control Language required. 
Basically, the calculation of the inviscid data tape 
is broken up into three steps as illustrated in the flow chart 
in Fig. 2. First the inviscid velocity distribution is gen- 
erated with either of the potential flow programs of Refs. 

3 or 4 . Next the surface coordinate system is generated by 
program BLOT. The inviscid velocity field is then interpolated 
onto the surface coordinate system, Fourier fitted and written 
on the boundary-layer code data tape in programs DERVAT and 

INVTAP. Program AOAT is used only with the axisymmetric HESS 
code and serves to combine the zero and ninety degree angle 
of attack solutions into the solution at the desired angle of 
attack. 
2. Description of Input Data 

Program HESS (Axisymmetric) 

A. Card Column Parameter Locations 


Card 1 - Header Card. 


-4- | 


ccl-60 
cc63-68 
Card 2 - 
ecl 
ec2 
cec3 
ec4 
ec5 
cec6 
cc7 
cec8-10 
eell 
ecl2 
ecl4 


ecis 
ecl6 
ccl7 
ecl8 


ecl9 


ec20 
ec21 
cec22 
cc23 
ec 28 
cce29,30 


Header 


Case 


Any alphanumeric run description. 


Case number. 


Control Flag Card 


NB 

NNU 
IAXI 
ICR@SS 
IQFF 
I@NLY 
IELPSE 
(blank) 
IPRTRB 
IPQTNL* 
IPTANV 


IVORT 

IQMITA 
LOMITC 
ISURFV 


IPRSCV 


ITALLV 
IEXCRS* 
IGENBC 
IRNGW 
IPNCH 
IUNIT 


Number of bodies (1<NB<9). 

Number of non-uniform flows (0<NNU<5). 
Axisymmetric flow flag. 

Crossflow flag. 

Off-body point input flag. 
Basic-data-only flag. 


Ellipse generator flag (see also Card 5). 


Perturbation velocities only. 
Solve potential matrix. 


Prescribed tangential velocity (for the last 
IPTANV bodies). 


Strip-ring vorticity flag. 
Omit axisymmetric uniform flow solution. 
Omit crossflow uniform flow solution. 


Surface vorticity (instead of sources) 
for the final ISURFV bodies. 


Prescribed values of the surface vortex 
strengths for the final ISURFV 
bodies will be input. 

All bodies are surface vorticity bodies. 

Extra crossflow. 

Generated boundary conditions. 

Ring wing option. 

Punched output. 


Unit number for input coordinates 
(default = 05). 


*Available if and only if N@NEWF = 1, ISIGF = 1 and IGEQMF = 1. 
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ec31 
ec32 
ec 33 


ec 34 


Card 3 - 
ecl-10 
ecll1-20 


cec21-80 
Card 4 - 
eel 


cec2 


ees 


ecé 


ec5 


cc6-10 


ecll-20 
cc21-30 
cc31-40 


cc41-50 


cc51-60 


IVIJ Matrix print flag. 


ICQEF Matrix-assembly coefficient print flag. 

IPRINT Very detailed matrix construction print 
flag. 

IRAKF Automatic rake generation flag (see 


also Cards 8 & 9) 
Chord/Mach number card. 
CH@RD Reference chord length (default = 1.0). 


XMACH Mach number for Goethert correction 
(0.0 implies incompressible). 


(blank) 
Body Control Card 1 of 2. 


IGEQMF 0 = curved elements; 1 = flat elements. 
ISIGF 0 = parabolic o; 1 = linear o; 
2 = constant o (on each element). 
ICURVN 0 = internally calculated element 
curvatures; 1 = input curvature 
(see card 7). 
NQNEWF 0 = use the newest formulas; 1 = use 


the old formulas (implies flat 
elements and constant ©). 


IFQRMT Input format flag (see Card Set 6). 


NN Number of defining endpoints for this 
body. 

XMULT x-multiplier value (default = 1.0). 

YMULT y-multiplier value (default = 1.0). 

THETA Coordinate rotation value (degrees, 


measured about -Z-axis. 


ADDX x-increment (to be applied to all the 
input coordinates for this body). 


ADDY y-increment (to be applied to all the 
input coordinates for this body). 
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Card 5 - Body Control Card 2 of 2 


ecl-10 IBDN "Body" number (sequential for bodies, 
zero for off-body points) 


ecll-20 ISUBKS Subcase flag 


ee21-30 NLF Non-lifting flag (for combination cases, 
only) 


ec31-40 A Semi-major axis for ellipse cases } If 

IELPSE 
ec41-50 B Semi-minor axis for ellipse cases } # 0 
Card (Set) 6 - Body Definition Cards 


IF IF@RMT=0: X-coordinates (6F10.5), then 
Y-coordinates (6F10.5) 


IF IFQ@RMT=1: X, Y coordinates (2F10.5) (i.e., one 
"point-set" per card). 


IF IF@RMT=2: X, Y coordinates (F10.5,10X,F10.5) 
(i.e., one "point set" per card). 


Card (Set) 7 - Input curvature values (needed only if 
ICURVN # 0) 


(€F10.5) CURV(I),I=1,NN-1 The curvature values for the 
NN-1 elements which constitute this 


Card 8 - Rake Number Card (needed only if IRAKF # 0) 


ccl-10 NRAKES The number of "automatically" generated 
mass-flow rakes 


Card 9 - Rake Definition Card (needed only if IRAKF # 0) 
eccel-10 Xi Coordinates of "start'' of the rake 
eeli-2z0 “YI 
cectoad #2 Coordinates of the "end" of the rake 
ec31-40 Y2 

ec41-45 N Number of intervals to be used in the 


rake (note, 4<N<200, and N must be 
an even integer) 


afm 


| Card 10 - Non-Uniform Flow Control Card (needed only if 
NNU # 0) 
ecl-10 NUN Flow identification number 
| 
F ecll-20 MSF 0 = axisymmetric onset flow; 1 = cross- 
| flow onset flow, 2 = both O and l 
ec21-30 TYPE +1.0 = velocity will be input in x,y 
component form; 0 = velocities will 
be input in normal, tangential form; 
-l = automatic generation of flow 
due to rotation about the Z-axis 
: (for crossflow, only) 
cc31-40 FG Flow generator constant 
Card (Set) 11 - Non-Uniform Flow Velocities (needed only 
if NNU # 0) 
(6F10.5) VX(I) or VN(I), I = 1, total number of 


control points 


i 

| 

(6F10.5) VY(I) or VT(I), I = 1, total number of 
control points 

| 


Card (Set) 12 - Specified Tangential Flow Velocities 
needed only if IPTANV # 0 


(6F10.5) TG(I), I = 1, total number of control 
points on the last IPTANV bodies 


B. Discussion of Parameters 

In all cases, the integer values (indicated by the standard 
FORTRAN IV naming convention) must be right adjusted within 
the specified input fields. For floating-point type input, 


the input decimal point will override the quoted FQ@RMAT specifi- 


cations. 


Card 1. 
Card 2. 


NB 


NNU 


TAXI 


ICROSS 
IOFF 


IONLY 


IELPSE 


IPRTRB 


IPQTNL 


Self-explanatory. 

Although most of these flags are unchanged from the 
original program, a brief description of each one 
will be given here for completeness. 

The total number of separate "bodies" (1<NB<8). 

The total number of (user specified) input non-uniform 
flows (0<NNU<5). 

Axisymmetric flow flag (0: no axisymmetric flow, 

1: axisymmetric flow) 

Crossflow flag (0: no crossflow, 1: crossflow) 
Off-body point input flag (0: no off-body points 
will be input, 1: off-body points will be input) 
Basic-data-only flag (0: full execution, 1: basic 
data only) 


Ellipse generator flag (see also card 5)(0: standard 


body defining coordinates will be input, 1: coordinates 


for the ellipse whose properties are specified on card 
5 are to be created automatically by the program). 
Perturbation velocities only (0: standard total net 
velocities including the onset flow will be printed, 
1: perturbation velocities only will be printed). 
Solve the potential matrix (0: standard solution for 
the surface velocities, 1: solution for the surface 
velocity potential). Since the velocity potential 
problem for the higher order elements has not been 


coded in the present version, the solution for the 
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IPTANV 


IVORT 


IQMITA 


IQ@MITC 


ISURFV 


potential is available only for the case of flat 
elements and, in particular, through the usage of 
the "old" velocity formulas. These "old" formulas 
have been retained as a subsection of this new, 
higher order program and may be "reached" by setting 
the following quantities (on Card 4): IGEQMF = 1, 
ISIGF = 2, and N@NEWF = 1; in such a case, the solu- 
tion for the potential is still available. 
Prescribed tangential velocity flag (0: standard zero 
normal velocity solution is to be obtained, 1: a 
user defined set of prescribed tangential velocities 
will be input on the last IPTANV bodies; obviously, 


IPTANV < NB). 


Strip ring vorticity onset flows for each of the 
bodies (0: no strip ring vorticity, 1: automatic 
generation of the strip ring vorticity onset flows). 
Omit the uniform axisymmetric flow solution (0: cal- 
culate the standard uniform axisymmetric flow solu- 
tion if IAXI # 0, 1: omit this uniform axisymmetric 
flow solution even if IAXI # 0). 

Omit the uniform crossflow solution (0: calculate 


the uniform crossflow solution if ICROSS # 0, 1: 


omit the uniform crossflow solution even if ICR@SS # 0). 


Surface vorticity for axisymmetric flow (in place of 


surface sources) flag (0: use the standard source 


elQ= 


| 
| 
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distribution on all elements, # 0: use vorticity 
instead of a source distribution as the singularity 
on the elements of the last ISURFV bodies). Note 
that if ISURFV # 0, then both IV@RT and IPTANV must 
be # 0, and then in particular, IPTANV must = ISURFV). 
If we assume that IPTANV = ISURFV = k, and that Mis 
the total number of elements on the last k bodies, 
then the surface vorticity option causes the source 
induced velocity formulas to be used on the first N-M 
elements of each matrix row, and the vortex formulas 
to be used on the remaining M elements of each matrix 
row. The solution for the unknown strengths then 
proceeds as it would have if only IPTANV was nonzero. 
IPRSCV Prescribed vorticity flag for axisymmetric flow (0: 
prescribed vorticity values will not be input, 1: 
prescribed vorticity values (W;) for the last IPRSCV 
bodies will be input; note that if IPRSCV is # 0, 
then both IV@RT and IPTANV must # O, and then, in 
particular, IPTANV must = IPRSCV). If IPRSCV is 


# 0, then the vortex strengths of the affected 
elements are taken to be W,/4n. 

IALLV Total vorticity flag for axisymmetric flow (0: all 
elements are assumed to be source type elements 
(unless ISURFV # 0), 1: all elements are assumed to 
be vortex type elements). Note that if IALLV is # 


0, the use of IPTANV is optional, but not mandatory. 


-ll- 


IEXCRS 


IGENBC 


Extra crossflow flag (0: no "extra" crossflow, 1: 
generate an "extra" crossflow (i.e., having a po- 
tential which varies as the cosine of twice the 
circumferential angle)). Note that if IEXCRS # 0, 
then NQ@NEWF = 1, ISIGF = 2, and IGEOMF must = l, 
since the higher order formulas for this kind of 
velocity potential have not been included in this 
program. 

Generated boundary condition flag (0: do not generate 
any additional boundary conditions for the crossflow 
case; 1: generate the onset crossflow due to rotation 
about an axis normal to the axis of symmetry (see 
also the notes for Card 10)). 


Ring wing option (0: no ring wing option, 1: use the 


ring wing option (see MDC Report J0741/01, April 


1970 for further details)). 

Punched card output flag (0: no punched card output, 
1: punched card output). 

This card is usually left blank, which results in a 
default chord length of unity, and no Mach number 
corrections (i.e., incompressible results). Ifa 
non-zero value of the Mach number is input, the pro- 
gram uses the Goethert correction to account for com- 
pressibility. 

Usage of the IGEOMF and ISIGF flags permit the user 


to "turn-off" any or all of the higher order element 
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Card 5. 


curvature and/or varying source density terms. 

The default values are curved elements with parabol- 
ically varying source density. If IGEPMF = 0, the 
program will automatically calculate the local ele- 
ment curvature values by the procedure described in 
Section II.A, unless ICURVN is non-zero, in which 
case the user must supply these curvature values 
(see Card 7). The value for the N@NEWF flag is 
ordinarily left blank (or zero), even if a flat 
element, constant source density solution is re- 
quired. However, since certain of the original 
Douglas-Neumann Axisymmetric Potential Flow Program 
capabilities have not been made available in the 
higher order program (e.g., calculation of the po- 
tential (hence the added mass, etc.) the original 
formulas have been preserved within this version (see 
description of the IP@TNL and IEXCRS flags). The 
original capacilities can be obtained by setting 
IGE@MF = 1, ISIGF = 2, and N@NEWF = 1 for each of 
the input bedies, in which case the input instructions 
of Reference 5 apply. 

The other parameters on Card 4 are self-explanatory, 
with the understanding that the order of coordinate 
transformations are as listed on the input. 

The "bodies" are normally loaded prior to any off- 


body points (although the latter require the usage 
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of body control Cards 1 and 2, also). For this 
reason, the value of IBDN should be sequentially 
increasing beginning with unity. A non-zero value 
for ISUBKS means that the body definition cards 

for this body (Card Set 6), (say this is the ith 
body input under this header card) will not be in- 
cluded, but that the program is to use the ith set 
of points that were input under the previous header 
card. Obviously, this capability is useful only 
with "stacked" input cases, 

The three input formats that are available are as 
shown. The value of IFMRMT determines which format 
is used. The default format is the "old" format: 
x-coordinates, followed by the y-coordinates. 

These values are to be input only if ICURVN # 0. 
This card is needed only if IRAKF # 0. Note that 


1 < NRAKES < 20. 


These are the rake definition cards, which define 

the "start" and "end"’ of each mass flow rake. The 

sign convention that is employed is as follows: 

a positive mass flux means that the flow is from left 
to right to an observer traversing the rake from 

point (Xl, Yl) to point (X2, Y2). For example, 

for an inlet without a bullet, point 1 is typically 

at y = 0.0, and point 2 would have a y-coordinate 
located on the inlet wall. Note that the x-coordinates 


need not be the same, i.e., a tilted rake may be used, 
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Card 10. 


Card 11. 


if desired. Note also that one rake may "overlap" 
another if so desired, since they are each treated 
independently. It should be pointed out that a 
Simpson Rule integration is performed over the ''N"' 
intervals (and therefore N must be an even integer), 
with the program automatically generating N-1l 


Las 


“intermediate” points as “off-body points" for each 
input rake (program limit is 500 total off-body 
points). The velocity values at the first and last 
rake points are obtained by linear extrapolation of 
the two nearest values. In this way, the difficulty 
associated with calculating induced velocities at 
off-body points which lie very near to, or on, the 
surface is avoided. Typically, values of N between 
10 and 20 appear to be satisfactory for most cases. 
Self-explanatory, except that FG, if entered is used 


in the following way to generate the rotation onset 


flow: 


. + si ft 
: COS oO; cP n as 


X. 


Self-explanatory. 


Program AOAT 


Card 1. Format (313,F10. 6) 
ecl-3 IPRT = 1, writes the output of program HESS 

on Unit 6; = 0, does not write 
cec4-6 N, Number of elements used in the program HESS 
cc7-9 KPL, Number of meriodonial planes required 
ecl0-19 ALPD, Angle of attack in degrees 


Program HFSS (3-D) 


Card 1. 
cc CONTENTS NOTES 
1+60 Title Header information 
61 IFLAG = "0" for input body coordinates (this 
: is the usual case) 
! = "1" for generated ellipsoids 
| = "2" for generated bodies with 
i elliptical cross section 
| 62 LIST = "0" for full execution 
| = "1" for basic data listing only 
| 63 KX = "F'" or blank for no generated uniform 
| x-flow 
= "T" for generated uniform x-flow 
(calculate x-flow matrices) 
64 KY = "F" or blank for no generated 


uniform y-flow 


= "T" for generated uniform y-flow 


(calculate y-flow matrices) 


=16- 


65 


66 


67 


63 


70+71 


74 


CONTENTS 


KZ 


ISIG 


IPRS 


MPR 


METH@D 


NNON 
NSYM 
NQOFF 


KMACH 


NOTES 


= "F" or blank for no generated 
uniform z-flow 

= "T" for generated uniform z-flow 
(calculate z-flow matrices) 

= "0" for no input guesses for source 
densities 

= "1" for input guesses for source 
densities (rarely used) 

= "0" for no mid-iteration print 

= "1" for print of solution iterations 
= "0" for no intermediate matrix print 
= "1" for be print 

a "2" fou AL. print 


ij 


= oO Vi; and Ay; print (Never use 


for large case, or output is excessive) 


= "0" for automatic solution selection 


"1" for direct matrix solution 

= "2" for modified Seidel iterative 
solution 

Number of input onset flows 

Number of symmetry planes 

= "0" for no off-body points 


= "1" for off-body points 


"O'' for no Mach number correction 


"1" for input Mach number 


eb 


ce CONTENTS NOTES 


75-76 == 


77-80 KASE Alpha-numeric case identification 


INPUT GENERATOR DATA 


Card 2. 
ce CONTENTS NOTES 
1. NLM1 Number of "latitudinal" element 
divisions 
6>10 MMIN Number of "longitudinal" element 
divisions 
11+20 B Y semi-axis of elliptical cross section 
21+30 C Z semi-axis of elliptical cross section 


(If ellipsoid is generated, X semi-axis is unity) 
NOTE: This card is input only when IFLAG # 0 
X-Z INPUT 
This input consists of the coordinate pairs (Xy; 21> 
Xo, Zo, - + - XypM1+1, Zu M1+D that define the profile in 
the x-z plane of the body of elliptical cross-section to be 
generated by the program. These coordinates are read with 


an 8F10.0 format. 


ON-BODY COORDINATES 


General (non-elliptical) bodies are input by giving the 
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coordinates of a large number of points on the body. The 
three dimensional coordinates defining the body surface(s) 
are input 2 points-to-a-card. Each point has associated 
with it a quantity called its "status" that has the follow- 
ing meaning: 

VALUE MEANING 
"0" or blank This point is on the same n-line as 


the last point 


1 This point starts a new n-line 
2 This point starts a new section 
3 This is the last point of the input 


NOTE: These cards are only inout if IFLAG=0 
OFF-BODY POINT COORDINATES 
The three-dimensional coordinates defining the off-body 
points are input with the identical format of the on-body 
coordinates. The only status flag used, however, is that of 
"3" for the last off-body point. 


NOTE: These cards are input only if NQ@FF=1. 


FLOW FLAGS 
Each input onset flow is preceded by a card that describes 
the flow. All x-flows must precede all y-flows which precede 
all z-flows. (An "x-flow'’ is any flow that is solved using 
the matrix appropriate for a uniform onset flow parallel to 


the x-axis.) 
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ee MEANING 

i) = "l' if this is an x-flow 
a "2" af thie is a ¥-flow 
= "3" 42 thie te a z-flow 

2 = "0" if this is a non-uniform onset 
flow, i.e., values of the onset flow 
velocities will be specified for each 
element 
= "lL" if this is 4 whiform onset flow, 
the value of which is only specified 
once 

a = "0" if the onset flow velocity is 
expressed in component form 
= "1" if the onset flow is expressed 
in normal velocity values 

NOTES: The onset flow velocities and sigma guesses are read 

with a "6F10.0" format. 


This input is only required if NNQN = 0. 


GENERAL PROCEDURES 
All coordinates are input with the following format 
Z (3Fi0.0, Ii) 
Except where noted elsewhere, zeroes or blanks may be inter- 


changed on input. 
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Program BLOT 


Card Name Fortran 


Col. 1-5 Name Format Column Description 
TITLE TITLE 9A8 6-77 Problem Identification - 


this will appear as a 
heading in the output 


LIMTS JMAX I5 6-10 Number of points across 
boundary layer 
KMAX 15 11-15 Number of points around 
body (in phi-direction) 
LITER E5 16-20 Number of iterations 
for initial profile 
IPRT 15 21-25 Every how many steps 
printing will occur along 
the body 
KPRT I5 26-30 Every how many steps 


printing will occur 
around the body 


PSIMAX E10.0 31-40 Distance along body at 
which program will term- 
inate | 
| 
UINF E10.0 41-50 Freestream velocity -V_ | 
DW DW(K) E10.0 6-15 Aw's around the body, 
K = 1, KMAX - l. 
16-25 
etc. 
DN DN(J) E10.0 6-15 An's across the boundary 
layer. 
J = 1, JMAX - 1. 
16-25 
SIZE A E10.0 6-15 a - length of axis of 


triaxial ellipsoid 
(A>106 gives paraboloids) 
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Card Name Fortran 


Col. 1-5 Name Format Column Description 
RB E10.0 16-25 R - radius of curva- 


ture of nose of body 
in the plane y'= 0 


RC E10.0 26-35 R, - radius of curva- 
ture of nose of body 
in the plane z' = 0 

EPS E10.0 36-45 Convergence criterion 


used in calculating the 


body coordinates 


ALP E10.0 46-55 a - angle of attack (deg) 


CAPTH E10.0 56-65 © - indicates which 


finite-difference scheme 


is being used. 


Z 


| 1 = Kruase scheme 
scheme 


CRI BLO. 66-75 


| 1 Implicit scheme 


1/2 Crank-Nicolson 


scheme 
Programs DERVAT and INVTAP 
These programs obtain the necessary data sets internally and 


they do not need any external input. 


3. Description of Program Output 
Program HESS(Axisymmetric) 
Unit 6 

The main output of the program, which consists of 
calculated surface velocities and/or pressures has not been 
changed. Its format is that of Ref.2. The surface coordinates 


with which the calculated velocities are associated are the 


a22= 


Crank-Nicolson 


transformed coordinates, i.e., the original input coordinates 
altered by any specified translations, multiplications, or 
rotations. Similarly, the format of the calculated velocities 
at the off-body points is similar to that of Pef.5.The only 
change is that output at rake points (q.v.) is included. 

The initial output, which details the surface geometry, 
has been changed considerably. The header identifies what 
kind of solution has been computed: (1) curved or flat 
elements; (2) constant, linear, or parabolic source density 


element curvatures. The header also specifies whether "new 

or "old" velocity formulas have been used, the latter of which 
apply only to the flat-element constant-source case. In the 
body of the output, the first two columns are the untransformed 
(input) coordinates, and the third and fourth columns are the 
transformed coordinates, which are the endpoints of the sur- 
face elements. The fifth and sixth columns are the coordinates 
of the control points of the elements. Element arc lengths 

are given in column 7 and a running total of arc length is 
displayed in column 8. Column 9 lists the differences between 
average slopes of successive elements, i.e., the differences 

in the slopes of the two straight lines, through the respective 
element endpoints. The actual slope discontinuity at an end- 
point between two parabolic elements is normally much smaller 
than the difference between the average slopes of the two 
elements. Moreover, this discontinuity approaches zero if the 
body contour approaches a parabola. The final (tenth) column 


lists elements curvatures. 


A new output is entitled Automatic Rake Calculation. 
For each rake the two input points that bound the rake and 
the input identification are output together with three calcu- 
lated quantities. The first is the surface area of the cone 
frustum defined by the two input rake points. The third is 
the total flux of fluid that crosses this area per unit time. 
The second quantity, average velocity, is the ratio of the 
third and the first quantities. There is a rake output for 
all axisymmetric flows, both uniform and non-uniform, but 
there is no rake output for cross-flow, because its circum- 
ferential variation guarantees zero flux. 

With regard to surface vorticity solutions, there are 
only two possibilities. Either there are no vorticity solu- 
tions or there is one solution for each body, which corresponds 
to a unit vorticity strength on that body and zero vorticity 
strength on all other bodies. The order of these solutions 
is the same as the order in which the bodies are input. Thus, 
in particular, an inlet with centerbody has two vorticity 
solutions - one with vorticity on the inlet and one with vor- 
ticity on the centerbody. The second of these is not meaningful 
and should be discarded. (It may, of course, occur first on 


the output.) 


Units l- 4, These are scratch disks, which are used 
7-15 to write and read internally. 
Unit 16 The coordinates x, y and the velocities 


Ty, Ty. Ty are written in order one 


after the other for all axial stations. 


atin 


Program AOAT 


Unit 6 


Unit 10 


Unte 22 


Unit 30 


Program BLOT 


Unit 6 


For each axial station (x), the pres- 
sure coefficient (CP), the velocities 
UX, UY, UZ are written. These velocities 
are in body axes. Finally, the stagna- 
tion point location XO is printed. 

For each axial station, two records are 
written unformated. The first record 
contains X and KPL. The second record 
contains the meridonial plane angle, 
UX, UY, UZ and radius R for all planes. 
This is not used. 


Formated image of Unit 10. 


Initially all input data is printed out. 
This includes several quantities which 
are computed strictly from the input. 
Next, initial profile calculations, in- 
cluding all iterations are printed. 
Finally, the solution at each desired 
location along and around the body is 
displayed. At each point where printing 
occurs, the following quantities appear: 

xX 

Y x,y,Z coordinates 


Z 
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r, radial distance to body surface 

s, distance along surface in &-direction 
ALPHA a = w,/u, 
ALBET2 = w,/u, 
ALBET3 = w./W 


XIKA parameter €& Ky /o 


UE u 
inviscid edge velocities in &€ and w 
WE Wo directions 


BETA1 By 
BETA2 Bo 

inviscid velocity gradient parameters 
BETA3 By 
BETA4 


DFDN 


Not calculated 
DGDN 


HPSI 


Metric coefficients 


HOMEGA 
KAPPAS 

Curvatures of the coordinate lines 
KAPPAW 


TOTAL SHEAR 
STRESS 


TOTAL SKIN Not calculated 
FRICTION 


ANGLE 


Unit 8 Formated image of Unit 25 
Unit 11 The diagnostic printout from the 
double interpolation subroutine SLOALL. 


Useful for debugging only. 
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Unit 25 


Program DERVAT 


Unit 6 


Unit 26 


Unit 29 


Unit 30 


For each (PSI) the longitudinal surface 
coordinate, the metric coefficients 
(HPSI, HW), R, S, UPSI, UW are written 


for each meridonial plane, unformated. 


The angle of attack in radians and the 
number of meridonial planes are written. 
Next the image of Unit 25 of program BLOT 
has been written. 

The angle of attack in radians and KPL 
are written. Next, the station No. I, 
XS, HPSI, HW, R, S, UPSI, UW, X, and 

the derivatives of these D(UPSI), D(UW), 
D(R), the second derivative D2(R); 
D(HPSI), D(HW), D(X), are written in 
order for each meridonial plane. These 
are written unformated. 

Formated image of part of Unit 26. The 
data written is I, XS, S, X, D(X), 

UPSI, D(UPSI), UW, and D(UW) for each 
plane. 


Formated image of part of Unit 26. First 


the angle of attack in radians and KPL 


are written. Next I, XS, HPSI, D(HPSI), 
HW, D(HW), R, D(R) and D2(R) are written 


for each plane. 
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Program INVTAP 


Unit 6 Frrors and messages are written on 
this unit. 
Unit 10 For each station the Fourier coefficients 


are written unformated. This is the data 
set required by ICBL3D. 

Unit 30 Formated image of Unit 10. Also at each 
station the data written on Unit 26 of 
program DERVAT is rewritten. 

JCL 

The programs described above require large memory, which 
made it necessary to execute them one after the other in a multi- 
step job. Also the programs were compiled, linkedited and stored 
as members of a partition data set, named A30303. MARK1. SUBLIB. 

This has the advantages of avoiding handling large Fortran source 
decks and compiling everytime the job is executed. The input 
cards for each program have to be inserted after the //FTO5FO01 DD* 
card, except for the program BLOT for which they have to be 

after the //FTO4F001 DD* card. Finally, the required inviscid 
data has been written on an on line disk (CHEKOV) as a partition 


data set named A5SO5SF3.INVTAP. The relevant JCL has been listed 


below. 

1 Job Card 

2 //JOBLIB DD DSN=A30303 ,MARK1 ,SUBLIB, DISP=(SHR,KEEP) , 
//VOL=SER=USERPK , UNIT=SYSDA 

3 //STEP1 EXEC PGM=HESS, TIME=1 , REGION=520K 

4 //FTO6FO01 DD DUMMY 
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//FTO1F001 DD 
=(TRK, (50,1)) 


//FTO2FO0O1 DD 
=(TRK, (50,1)) 


//FTO3F001 DD 
=(TRK, (50,1)) 


//FTO4FOO1 DD 
=(TRK, (50,1)) 


//FTO7FO0O1 DD 


//FTO8FOO1 DD 
=(TRK, (50,1)) 


//FTO9FOO1 DD 
= (TRK, (50,1)) 


//FTLOFOO1 DD 
=(TRK, (50,1)) 


//FT11F001 DD 
=(TRK, (50,1)) 


//FT12F001 DD 
=(TRK, (50,1)) 


//FT13F001 DD 
=(TRK, (50,1)) 


//FT15F001 DD 
=(TRK, (50,1)) 


//FT16FO0O1 DD 


=(TRK, (50,1)), 


//FTO5FOO1 DD 


UNIT=SYSDA , DISP= (NEW , DELETE) , SPACE 


UNIT=SYSDA , DISP=(NEW, DELETE) , SPACE 


UNIT=SYSDA, DISP= (NEW, DELETE) , SPACE 


UNIT=SYSDA , DISP=(NEW ,DFLETE) , SPACE 


DUMMY 
UNIT=SYSDA, DISP=(NEW, DELETE) , SPACE 


UNIT=SYSDA , DISP=(NEW, DELETE) , SPACE 


UNIT=SYSDA , DISP= (NEW, DELETE) , SPACE 


UNIT=SYSDA,DISP= (NEW, DELETE) , SPACE 


UNIT=SYSDA, DISP=(NEW, DELETE) , SPACE 


UNIT=SYSDA , DISP=(NEW, DELETE) , SPACE 


UNIT=SYSDA, DISP= (NEW, DELETE) , SPACE 


UNIT=SYSDA , DISP=(NEW, PASS) , SPACE 
DSN=&&T1 


* 


//STEP2 EXEC PGM=AOAT,COND=(0,NE) , TIME=(0. 30), 


REGION=150K 
//FTO6FOO1 DD 


//FT16F001 DD UNIT=SYSDA,DISP=(OLD, DELETE) , DSN=&&T1 


//FT22F001 DD 
//FT30F001 DD 
//FT1OFOO1 DD 


DUMMY 


DUMMY 
DUMMY 
UNIT=SYSDA, DISP= (NEW, PASS) , SPACE= 


(TRK, (50,1)) , DSN=&&T2 
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25 //FTOSFOO1 DD * 

26 //STEP3 EXEC PGM=BLOT ,COND=(0,NE) , TIME=5 ,REGION=350K 

27 //FTO6FO01 DD SYSOUT=A 

28 //FTO8FO01 DD SYSOUT=A 

29 //FT11F001 DD DUMMY 

30 //FT1OFOO1 DD UNIT=SYSDA,DISP=(OLD , DELETE) , DSN=&&T2 

31 //FT25F001 DD UNIT=SYSDA,DISP= (NEW, PASS) , SPACE= 
(TRK, (50,1)) , DSN=&&T3 

32 //FTO4FO01 DD * 

2 //STEP4& EXEC PGM=DERVAT, COND=(0,NE) , TIME=1 ,REGION=650K 

34 //FTO6FOO1 DD SYSOUT=A 

i //FT25FO01 DD UNIT=SYSDA,DISP=(OLD, DELETE) , DSN=&&T3 

36 //FT26FO01 DD UNIT=SYSDA, DISP=(NEW, PASS) , SPACE= 


(TRK, (50,1)) , DSN=&&T4 


37 //FT29F001 DD SYSOUT=A 

38 //FT30F001 DD SYSOUT=A 

39 //STEPS EXEC PGM=INVTAP ,COND=(0,NE) , TIME=1 , REGION=650K 
40 //FTO6FO01 DD SYSOUT=A 

41 //FT25F001 DD UNIT=SYSDA, DISP=(OLD, DELETE) , DSN=&&T4 

42 //FTLOFOO1 DD UNIT=SYSDA,DISP=(NEW,CATLG) , SPACE= 


(TRK, (50,1)), 
// DSN=A505F3.INVTAP , VOL=SER=CHEKOV 


43 //FT30F001 DD SYSOUT=A 


JCL for use with Fortran source decks 
In the JCL listed above the following changes have to 
be made: 


1. Card 2 has to be removed. 


~Siie 


ul 


2. In card 3 (PGM 


HESS) has to be replaced by (FORTGCG) 


19 ¢ = AOAT) has to be replaced by (FORTGCG) 

26 ( = BLOT) has to be replaced by (FORTGCG) 

33 ¢ = DERVAT) has to be replaced by (FORTGCG) 

39 ¢ = INVTAP) has to be replaced by (FORTGCG) 
3. After cards 3, 19, 26, 33, 39 insert the following 


card: 


//FORT.SYSIN DD* and then appropriate Fortran source 
deck. 


Listings of Job Control Language, sample data decks, 
sample output and program listings for the TAPGEN package are 


found in Appendices I through V. 
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SECTION II 


PROGRAM ICBL3D 


In this section program ICBL3D is described. The 
general structure of the program is briefly described and 
then the input data and output of the program are described 
in detail. Finally, the Job Control Language required to 
execute the program on an IBM 370/158 computer with the JES2 
operating system is listed. 

A flow chart of the program appears in Fig. 3 where 
the function of the various subroutines are indicated. As 
can be seen, the code requires a small amount of card input, 
but its' primary mode of data input is from a data tape desig- 
nated as logical unit 10. This data tape can be prepared 
using program TAPGEN described earlier in this report. 

The program proceeds in a forward marching fashion 
from the body stagnation point downstream and from the windward 
to leeward symmetry planes. The number of solution planes around 
the body in the w-direction is fixed by the user, while the step 
size in the €-direction is adjusted internally by the program 
as it proceeds downstream. The initial step size is set by the 
user and this step size is continually adjusted by the program 
based upon the number of iterations required to achieve a solu- 
tion at a previous station. Further, at the user's option, 
the edge value of the transformed normal coordinate 1 can be 
automatically increased or decreased to meet a predetermined 


criteria for profile decay at the boundary-layer outer edge. 


32% 
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There is a special feature of the program for full 
three-dimensional solutions. If the program cannot obtain a 
converged solution at some point between the windward and lee- 
ward symmetry planes it will drop that point (and all subse- 
quent points to the leeward plane) from the solution around 
the body at all subsequent streamwise stations. This feature 
allows the solution to proceed even after separation or other 
problems are encountered near the leeward plane of the body. 
1. Program Input 

As mentioned, input to the program is from two different 
sources: (1) cards or card image data and (2) data sets re- 
siding on tape or disk. The data sets residing on tape or 
disk are generated by program TAPGEN and basically consist of 
Fourier coefficients. The first record on the tape contains 
the angle of attack in radians and an integer variable giving 
the number of coefficients ineach Fourier series. All the 
records on the file after the first are identical in format. 
They consist of the value of & followed by the Fourier coef- 


Aicients for the series fits to r; U.,Wi,h., bh, Ss, dU/0e,; 
e Se gt ce) e 


aw /36, ar/ae, a°r/aey, ah,/3E, 2h,/9E, x, ax/3E in that order. 


The series are in terms of the angle > and their range is from 
the windward to leeward planes of symmetry for constant & 
All data on the tape are written with unformatted WRITE's. 
The data must appear on logical unit 10. 
The card or card image data to the program will now be 


described. 


2. Description of Input Data 
CARD 1. LABEL (20A4) 

LABEL is the title of the case and is a single sub- 
scripted array. LABEL is passed to the output routines and 


appears as the header to the program output. 


CARD 2. IE (49X, 13) 

IE is the number of points taken normal to the 
body. The program is dimensioned for a maximum of 101 points 
in the normal profile arrays. 101 is the recommended value; 
however, savings in execution time occur by decreasing IE. 


IE must be an odd number. 


CARD 3. INJCT (49X,13) 
INJCT is the subscript of the XSTA array (see below) 
giving the surface location at which mass injection begins. 


A zero value of INJCT is set to NSOLVE. 


CARD 4. KADETA (49X, 13) 

KADETA is an indicator for the adjustment of the 
transformed normal coordinate. If KADETA = 0, the maximum value 
of n is held constant. If KADETA = 1 the maximum n is adjusted 
when the velocity profiles fail to decay properly at the outer 
edge as prescribed by the input variable ADTEST. The value of 
n as adjusted up or down for windward symmetry plane only 


max 


runs. For full three dimensional runs Nive is only adjusted up, 


and can be adjusted at every point. A value of 1 is recommended. 


CARD 5. KEND2 (49X,13) 


KEND2 specifies the number of circumferential planes 
to be used in the solution, and thereby also sets the circum- 
ferential step size. The program can accept a maximum value 
of KEND2 = 61. A total of sixty-one planes would lead to 
excesSive computing time, however, and a value of KEND2 = 11 


is recommended. 


CARD 6. KONSET (49x, 13) 

KONSET is the subscript of the XSTA array giving the 
location of the onset of transition. At X = XSTA(KONSET) 
the variable LAMTRB is reset to 2 and transition to turbulence 


begins. A zero value is reset to NSOLVE. 


CARD 7. KPRT (49X,13) 

KPRT is a print control parameter which controls the 
printing of profiles in the ¢ direction. The program prints 
every KPRT'th profile around the body at a given value of & 


starting with the ¢ = 0° profile. 


CARD 8. KTRANS (49X, 1S) 
KTRANS is an indicator for the transition model. Bi 


KTRANS 


O transition to turbulence will be instantaneous, if 


KTRANS 


1 a smooth transition to turbulence will take place 


over the distance XBAR (see below). 


CARD 9. LAMTRB (49xX,13) 
LAMTRB indicates whether the flow is laminar or tur- 


bulent. LAMTRB = 1 indicates the solution begins with laminar 
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flow. LAMTRB = 2 means fully turbulent flow. When KONSET 


# O, LAMTRB = 1. 


CARD 10. LPRT (49X,13) 
LPRT is the print control parameter in the streamwise 
direction. The program prints solutions at every LPRT'th 


station in €&. 


CARD 11. NIT1 (49X,13) 

NIT1 is an iteration counter used to adjust the stream- 
wise step size. If the total number of iterations required to 
obtain the solution at a point is less than or equal to NIT1 


then the & step size is doubled. 


CARD 12. NIT2 (49X,13) 

NIT2 is an iteration counter used to adjust the stream- 
wise step size. If the total number of iterations required to 
obtain a solution at a point is greater than NIT1] and less than 


NIT2 the & step size is unchanged. 


CARD 13. NIT3 (49X,13) 

NIT3 is an iteration counter affecting — step size and 
the convergence of the solution. If the total number of itera- 
tions required for a solution at a particular point is greater 
than NIT3 the program halves the € step size and cuts back the 
value of & by the new step size. A solution at the smaller 
value of € is then tried. If this procedure fails three con- 
secutive times execution is terminated. 


NOTE: The & step size is adjusted only at the windward 
symmetry plane. 
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CARD 14. NOINJ (49X,13) 
NOINJ is the subscript of the XSTA array giving the 
surface location at which injection ends. NOINJ = 0 is reset 


to NOINJ = NSOLVE. 


CARD 15. SFC (49X,A3) 
SFC is a literal variable coded as either YES or NO 
to indicate whether surface curvature effects are to be in- 


cluded in the calculations. 


CARD 16. NSOLVE (49X, 13) 

NSOLVE is the number of variables in the XSTA array, 
and is therefore the subscript of the last XSTA value which 
indicates the end of the body. It is also the default value 


for INJCT, NOINJ, and KONSET. 


CARD 17. ADTEST (49X,E14.6) 

ADTEST is used in conjunction with KADETA. When KADETA 
is 1 ADTEST provides the convergence criteria for checking the 
streamwise velocity profile. When |F(IE) - F(IE-4)| or |G(IE) 
-G(IE-4)| is less than ADTEST/10 the maximum value of 1 is 
decreased by 10%. When it is greater than ADTEST the maximum 


value of nis increased by 10%. 


CARD 18. AKSTAR (49X,E14.6) 
AKSTAR is a numerical constant in the Van Driest inner 
eddy viscosity law (k* in Ref. 1). The recommended value is 


0.435. 


tJ < 


CARD 19. ALAMDA (49X ,F14.6) 
ALAMDA is a numerical constant in the outer eddy 
viscosity law used in the program (A in Ref. 1). The 


recommended value is 0.09. 


CARD 20. ASTAR (49X,E14.6) 
ASTAR is a numerical constant used in the damping term 
of Van Driest's inner eddy viscosity law in the program (A* in 


Ref. 1). The recommended value is 26.0. 


CARD 21. CWALL (49X,E14.6) 
CWALL is the injection rate for mass transfer cases 


where the rate is a constant. CWALL = Be 


CARD 22. CONV (49X%,E14.6) 

CONV is the solution convergence criterion. The de- 
pendent variable arrays of streamwise and cross flow velocities 
are all checked for convergence at all points. When the largest 
percentage difference between the current and previous iterations 
is less than or equal to CONV the solution is taken to be con- 


verged. 


CARD 23. DXMAX (49X,E14.6) 
DXMAX is the maximum step size, to be taken in 


the streamwise direction. 


CARD 24. DX1 (49X,E14.6) 
DX1 is the initial streamwise step size. Since 


this value will be adjusted internally it is not critical that 
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the user choose an accurate value. Usually a value of 0.01 


is a good initial DX. 


CARD 25. EDYLAW (49X,A3) 

EDYLAW specifies the inner eddy viscosity law to be 
used in turbulent cases. Two options are available to the 
user: (1) EDYLAW = VAN DRIEST and (2) EDYLAW = REICHARDT. 
The program picks up only the first 3 letters of each name. 
The user should note the comments on these two laws for mass 
transfer problems stated in Ref. 1. In general the REICHARDT 


law is recommended. 


CARD 26. ETAFAC (49X,E14.6) 


ETAFAC controls the normal grid spacing. A value of 


a i th nokta sl 


1.0 gives an equally spaced grid for the transformed normal 
coordinate. A value greater than 1.0 gives a finer grid at 
the wall than at the outer edge. A value of 1.04 is recom- 
mended with IE = 101 and ETAINF = 6.0. A value of 1.09 is 


recommended with IE = 101 and ETAINF = 100.0. 


CARD 27. ETAINF (49X,E14.6) 

ETAINF is the maximum value of yn. A value between 6.0 
and 10.0 is recommended for laminar flow. A value between 10.0 
and 100.0 is recommended for turbulent flow. This value can be 


adjusted internally by specifying KADETA = 1. 


CARD 28. PR (49X,E14.6) 


PR is the value of the laminar Prandtl number in the fluid. 
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CARD 29. RTW (49X,E14.6) 
RTW is the ratio of the wall temperature to stagna- 
tion temperature. It is used to calculate wall temperature 


when wall temperature is held constant. 


CARD 29. TFS (49X,E14.6) 
TFS is the free-stream static temperature in degrees 


Rankine. 


CARD 30. cP (49X,E14.6) 


CP is the specific heat of the fluid in £t7sec?°R. 


CARD 31. AMUINF (49X,E14.6) 


AMUINF is the coefficient of viscosity in slugs/ft-sec. 


CARD 32. PSTAG (49X ,E14.6) 


PSTAG is the free-stream total pressure in bf /ft?. 


CARD 33. PINF (49X,E14.6) 


PINF is the free-stream static pressure in lbf/ft?. 


CARD 34. XBAR (49X ,E14.6) 
XBAR is the relative length of the transition regime 
in turbulent cases. It is the ratio of the transition end 


point location to the transition onset distance. 


CARD 35. UFS (49X,E14.6) 


UFS is the free-stream velocity in ft/sec. 


CARD 36 to 36+NSOLVE XSTA(I), I = 1, NSOLVE (Fi2.6) 
XSTA is an important input array. It is an array of 


surface distances 


EO AEN A tA Tp COLD ERLE ELIE SI I et 


solutions calculated. The program will always obtain solutions 

at these points regardless of internal adjustments to the stream- 
wise step size. Both the value 0.0 and the end point of the body 
must be included in the array as well as any distances describing 


the beginning or end of injection and transition. 


CARD 36+NSOLVE+1 TO LAST CARD XTW(I), TWX(I), XCI(I), CIX(I) 
(4E12.6) 


TWX(I) is the wall temperature in degrees Rankine at 
XTW(I) which is a surface distance. 

CIX(I) is the injection rate ¥/8., at XCI(I) which is 
a surface distance. 

These arrays are dimensioned for a maximum of 500 values 
each. If none of these cards appear in the input deck the pro- 
gram will automatically assume constant wall temperature and 
injection rate values based on RTW and CWALL. This input allows 
the wall temperature distribution and injection rate distribu- 
tion to be read in versus their own surface distance tables. 

If both distributions are to be read in versus the same distance 
table, then either one of the two distance tables may be left 
blank. Another important feature is the fact that the distri- 
butions need not cover the same surface distance. For instance, 
the wall temperature distribution might be defined over the 
entire body while the injection rate distribution might only be 


defined over a short distance. 


3. Description of Output Data 
The output from program ICBL3D is in the form of 


printed output. The printed output is presented in two forms, 


me 


I 


station data and profile data. This section will define the 
output by listing the variable name as it appears in the output 
along with the definition of the variable. The program also 
prints miscellaneous messages which are described at the end 


of this section. 


Station Data 
Line l 


XO 


PHI 


Line 2 
XI 

DXI 
CWALL 
NIT 
Line 3 
HX 

HW 
Line 4 


TE 


distance from the stagnation point along 
the body axis non-dimensionalized by the 
reference length (one foot) 

radial distance from the stagnation point 
body axis to the body surface non- 
dimensionalized by the reference length 
meridional angle about the stagnation 
point body axis measured from the windward 


symmetry plane in degrees 


longitudinal surface coordinate & 
A&, streamwise step size 
local injection rate %,/.. 


number of iterations to obtain the solution 


metric coefficient h- evaluated at the body 


E 


metric coefficient ho evaluated at the body 


nondimensional edge static temperature, 


2. 
ToC p/ Yeo = T, 


UE nondimensional & component of velocity 


at the boundary-layer edge u./u,, =U 


e 
VE nondimensional w component of velocity 
at the boundary-layer edge ¥,/U,, = Re 
Line 5 
DUEDX dU ,/9& 
DVEDX oW/98 
DUEDW sU.,/ dw 
DVEDW oW/3u 
Line 6 
LOCAL EDGE REYNOLDS NUMBER = PU QE/U,, 
Line 7 
CFXINF Cc = 2 /oe fo streamwise skin friction 
aR G 00 “00 ae 
ee Wy coefficient 
CFXEDG Cr = 2. me streamwise skin friction 
be We coefficient based on edge 
conditions 
CFWINF Ce = 2. jena, transverse skin friction 
We Ww coefficient based on 
free-conditions 
CFWEDG Cc = 2 /0 o transverse skin friction 
af T oe Stas 
We Wo coefficient based on 
edge conditions. 
Line 8 
QW dy/P.U2 wall heat transfer coefficient 
CHIMAX maximum vorticity Reynolds number n2/v du/dn 
Line 9 


LONGITUDINAL SKIN FRICTION a in lbf/ft? 
x E 
DELTA*(X) Ses streamwise boundary-layer displacement 


thickness in feet. 


oh» 


THETA(X) 


Line 10 


TRANSVERSE SKIN FRICTION 


DELTA*( PHI) 


THETA(PHI ) 


Line 11 


WALL HEAT TRANSFER RATE 


DELTA(FT) 


Profile Data 


A) streamwise boundary-layer 


&? 


momentum thickness in feet. 


Tyo in lbf/ft? 


Se transverse boundary-layer 
displacement thickness in feet. 
ie transverse boundary-layer 


momentum thickness in feet. 


a. in Btu/ft~/sec 
6, the boundary-layer thickness in 


feet. 


One group of boundary layer profile data is printed by the 


program. Every other point is printed in the profile arrays. 


The following columns are identified: 


ETA n, the transformed normal coordinate. 


Y n, the physical normal distance nondimensionalized 


by the reference length. 


F h-u/hy QU, 
FN dF/3n 
G h w/h ow where W = * at the stagnation point and 


symmetry planes and W = U, elsewhere. 


GN dG/an 

- T/T, 

TN 9(T/T,)/an 

V transformed normal velocity profile. 


EPLUS er the eddy viscosity coefficient 
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Aa ase 


Miscellaneous Messages 

The program has a few internal messages which are 
written to indicate problems with the solution, or coordinate 
adjustments. A message is printed by subroutine ADDETA when- 
ever n, is adjusted up or down. The direction of adjustment 
is given along with X, the old n, and the new n,. 

A message is printed by subroutine CHANGX indicating the 
beginning of transition or mass transfer. Included in the 
messages are the values of X and the particular integer counter 
involved. A similar message is also printed by CHANGX when mass 
transfer ends. 

Whenever the program fails to obtain a converged solu- 
tion within NIT3 iterations, a message is printed by subroutine 
CONTRL to that effect which includes the values of the trans- 
verse and streamwise solution counters and NIT. If this should 
occur turee consecutive times, a message will be printed indi- 
cating that execution is terminating. 

If a particular boundary-layer problem drops all of its 
circumferential solution planes due to convergence problems, a 
message will be printed by CONTRL indicating that execution is 
terminating. 

A normal termination of the program is indicated by 


the message "THE END" printed out after the last station results. 


3. Sample Input/Output 
Listings of Job Control Language, sample input decks, 
sample outputs and format sheets for the input data cards and 


the program listings are found in Appendices VI through X. 
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Figure 1. BODY COORDINATE SYSTEMS FOR INVISCID 
AND BOUNDARY-LAYER CALCULATIONS 
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Axisymmetric HESS 
a = 90° 


a =90° HESS solutions 
for solution at 
correct «a 


BLOT 
0° and 


0° and 


DERVAT 
Fit data in &-direction 
and calculate pressure 

gradients 


INVTAP 
Fourier fit geometry 
and inviscid data and 
write on tape 


ICBL3D 
Predict body boundary- 
layer 


BLOCK DIAGRAM FOR ICBL3D AXISYMMETRIC BODIES 


ahJ= 


Generate coordinate 
system and metric 
coefficients 


3-D HESS 
Calculate inviscid 


BLOT 


Figure 2b. 


DERVAT 
Fit data in &-direction 
and calculate pressure 
gradients 


INVTAP 
Fourier fit geometry 
and inviscid data and 
write on tape 


ICBL3D 
redict body boundary- 
layer 


BLOCK DIAGRAM FOR ICBL3D NONAXISYMMETRIC BODIES 
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Generate coordinate 
system and metric 
coefficients 


Figure 3. 


READ INPUT CARD DATA 
(INPUT) 
INITIALIZE VARIABLES 
(INIT) 


CALCULATE Re 
(AERO) 


CALCULATE BOUNDARY LAYER 
(CONTRL) 

BEGIN STREAMWISE 
ITERATION LOOP 
BEGIN CIRCUMFERENTIAL 
ITERATION LOOP 


OBTAIN EDGE VALUES AND 
METRIC COEFFICIENTS 
(EGPROP , BLUNT2 , HOFN, GEOCOF , GGCOF) 


SOLVE STREAMWISE MOMENTUM 
EQN 
(XMOM, ABCDE, SOLVE) 


CALCULATE TRANSFORMED NORMAL 
VELOCITY 
(VCALC) 


SOLVE CROSS FLOW MOMENTUM 


EQN 
SUBROUTINE (PHIMOM, ABCDE, SOLVE) 
FLOW CHART OF PROGRAM ICBL3D 


«&9- 


COMPUTE TRANSFORMED 
NORMAL VELOCITY 
(VCALC) 


DROP CURRENT 
PLANE FROM 
SUBSEQUENT SOLNS 


YES (K=1) VES (K>1 


YES 


TERMINATE 


EASCURECH CONVERGENCE TEST 


PASSED 


FAILED 


CALCULATE B.L. PROPERTIES 
AND PRINT SOLUTION 
(PROPTY, OUT2) 


RESET PROFILES AND COUNTERS 
FOR SOLUTION AT NEXT POINT 


LAST CIRCUMFERENTIAL NO G) 
PLANE? 


YES 


CALCULATE NEW & 


TEST FOR BEGINNING OR END 
OF MASS TRANSFER AND 
TRANS ITION/ TURBULENCE 


@) ae LAST £2 
YES 


END 


Figure 3 continued: FLOW CHART OF PROGRAM ICBL3D 


( ) = Subroutine Names 
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APPENDIX I 


Job Control Language for Program TAPGEN 


OS/17/7TT 


DATE 


520K,PA2M.F IRT="NOSOURCE® 


4, REGION 


foe 


FORTRAN SOURCE DECK 


PROGRAM HESS 


//GOeFTOSFOOL OD * 


HESS 


PROGRAM 


INPUT DATA DECK FOR 


eeeree ere © © & Oe 
lelelelelelolojojeloleo! |) 
LALA LAU LALA LALA oe 
ee ie) 
ceeereeerteeeeeyy 
<~x¥¥YY¥YYtYYYyYY~ 

[odie Me de Go Aen Ao a Ae Aa ae 

cell eee ane eel eee coed el cool coll el eel 
ee we ey 

Hn nh ww ee 

Us a dd 
OVUOUVULOULUL UDO II 
adqdaqddaqdaqddw 
aaaaaagaanadayw 
NANANNNANNNNANY 
eeeeeereee & of 
eae ien hee tet OA) 
WWW Lae 
i 
WWW Wi 

as es De Ds ee ee ee ee ee) 
WWW a a le 
OCO000CRe0O00a 
eeereerrteeree ee 
ESET STSETSTETIET 

WO LOU LW Ws LW Ll 
222222222222 
ee we ee ewe 

Ce 
aaqaaagaagqagagagaang 
NNNNNMNMNNNNNWY 
tk ee ad aed eet ed nd pee et beg De 
QOAQCCHOOCA00A0G 
eorerererere en ene 
ddaqaqddddqdaa 
AOAODQ0OCRhRA0684 

A NNNNANNMWMNMNANNMN 
ral >>> >> >>> >>> > 
NDNANANNNMMNNNIWY 
S> hn wn 
[0 > al coal cell eel ol eel col aed el al lo 
(Va > We el heel a) 
D222 222222222 
AGDDID> D255 55 


20000000900000 
BOotoO0CHaCOoCaA vo2 


L,REGION=300K 


=(0,NE),TIME= 


ORTGCG,COND 


Nd el a ed to od ted OC /) 
lolelslololelolelelolelelolome*) a 
O000000090900T099 ~n 
UUUUUUUUUUNUEn Ow ® 
OMAN POROANMUNO Ae 
lelelelelelelolel ltt femme. 4 
segue ened ane] eel and eoelh au ened aed anal ceed saad eee 
[oS ES ey ea a Ye 
SAAR RRR RRS ESS 
SSwe~nfe~n few n~newnn~»w~nBwn 


FORTRAN SOURCE DECK 


PROGRAM AOAT 


//GO.FTOSFOOL OD * 


AOAT 


PROGRAM 


DATA DECK FOR 


INPUT 


=GETL 
=500K 


10,REGION 


N 
ws 
us 
" 
z 
wn 
a 
o 
~ 
~ 
~ 
o 
jo) 
w 
e 
x 
[<4 
- 
~ 
" 
us 
oO 
a 
oO. 
Y 
o 
wn 
n" 
<= 
a 
= 
pee 
2 
~ 
" 
a 
WY 
_ 
=) 
co 
<a 
a 
n 
> 
Za) 


SYSDA,DISP=(OLD,DELETE) ,OSN 


FORTGCG yCOND=(0,NE), TIME 


~n NSN ENS 
~nN S 


FORTRAN SOURCE DECK 


PROGRAM BLOT 


//GQO.FTOSFOOL DD DUMMY 


BLOT 


PROGRAM 


INPUT DATA DECK FOR 
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| 
| 
I 


O5/17/T7 


DATE 


) 
) 
(TRK (50,19), 


R 
R 
) 
CE 


aawaayt 


E 

E 

L 

L 

E 
4,REGION=500K 


D 
D 
V 
D 
) 
T 


ANNCKANL 
rr nnm tnd 
SEEZer ee 
TEL RH OH y 
DDDZZ022NY 
S0o0D9D8D750* 


QaC000 C0 
Oacaagq 00 


Heat td 
ooo000 0° 
lolelelelo mmole} 
WU Wi 
ODDAINA OW 
COaaHt aN 
RRR ee ee 
ST oe AW 
.SSax SARA SENN 
sSNA SSS SR SSNS 


4 EXEC FORTGCG,TIME 


TEP 
ORT.SYSIN DD * 


FTO4FO001 DO 


FORTRAN SOURCE DECK 


PROGRAM DERVAT 


=EET4 


1)),DSN 


Pe 


FOl 
=(T 
=300K 


S 
P 
4,REGION 


O E 
N S 
(O,NE),TIME= 


A,D 
A,D 
OND= 


FORTRAN SOURCE DECK 


PROGRAM INVTAP 


E=(TRKy(50),RLSE), 
USERPK,y 
=EET4 


p 
E 
TE) »OSN 


KIN INT 
ZWuKDOoY I 
FRANNWHE 
.DdniaeernD 
OOF seo 

Nm ZF omy 


sSARNNN EN 
sSSAnxxnn SNA 


=Shn 


APPENDIX II 


TAPGEN Sample Input 


s2°0 o°t 


1 
SiN3nW3T3 O21 930 OT =VHdIV¥ i 
SS3H wvHd0Nd 


L4L/12/40 310 


NNN 
DOCOMMMA MAM 
ee eee esee 


NANN 
oooanmmman 
2 0 2S 4-8 (6 6 


NAN 
soommamm 


ANGLE OF ATTACK 


SNNNN 

Qoo0commmann 
oe 

0 


8LOT 


-0000001 +10.0 


“oO0COmMmMmAmAN 
oer eeetaeve 


iol le lelelelslslalalia 
mH Oe eee eee ee 


“9 
ZNAOSOCOMMMMMOO 
maIme tee eee eee 


O4/21/T7 


DATE 


APPENDIX III 


TAPGEN Sample Output 


LUTION 


en 
DATA 


FLW FIELL 


NOITIONS FAR ANUNTARY LAYER 
INVISCID 


FROM Tht 


co 
TAKEN 


INVISCIO EDGE 


PA 


INVISCIO DATA 


IN THE 


NUMBER OF PLANES 


ALPHA= 


| 


ere 


LOw 


F 


INVISCIY 


OLUTION FOR THE 


cS 


UNTFORM FLOW STARTING 


WALL DATA AT STATION 


tt 


= 


= | 


UPS! 


Hw 


HPST 


K 


Lr oe tie pe eer ek 
oer eevee toes eoe 
Sa we SIO Ch Prete ICH! 


a had EAIOBoSasct = ~~ Wes FO ry hie: 
oe e¢ee 8 € @ € 8 C.8.8 € CE OSS B18 6 EO @ 
JWL OU VOCE UF UO IO CCU VL = po Se 


. 
OOCCOGCLVOGC GC 


VOovSEeCGcooccur VS 
eoveoeeceeree ee eevee 


eccooscecooocscocoococococoe 
ee 
CETOOCOOCCE OSOOC COCCOCUSOS 


OCOT MOQODOCOCOOGCOCOC OOOCTC OOCOGEL 
eeeeeeaeteveseveteeeeeeneeeteeneeae 


ECCECOTOCDOCOCC EC SOSC CO VOGECOVSe 


Et tt tt tt tt i et tt 
BOCOCCOSCCGOGCOSTCOOOOOTC GCOO0G 


} cS 5° 
SECOSCOSOCCOCoOCSCo CoO SC SCOCSO 
tt tt te et tt tt et tt tt Et et ed 
eee ee eee een e eee eeeeeeeeeee 


COOCTOCOTCOOCOOOCSOOSCOOOCOC COO 


ANOS HOME CORAM INOFDOORNM IMHO 
A Rt mt MMA 
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O (Hw) x) 


OR) H2(R) UCHeST) 


OtUw) 


O(uPSI) 


K 


Et et tt tm et ie 


IDE OL ODS OWIDIOGSO 
itn ea ice i tat 
J IO 23 23- 
JSHRISITOCI Sane mer a ay 
VO NAN NANNN 
‘ Vv oItV 

DT 6 
evens 
Aor 
see 
COUs 

cs 


elt 


Souse ¢ UUIGUCO SEC Se 
PIIITSIRGG sus 

3D 264. 00e C2 
tttersrsts It 
~~ SRA NR re 
Smee ae om 
eit cCscgle ce 
oerteeeeeee ee 
IIE SICIIIG om 


Ol ell el cll nd ad el ell el eet 
HS UVCCOUDGCEOUCCeED. 


Jot RR aH eta NOU 
MNLF HA wT VA DADMNOADM: 

c SS OMGOCOOU GO. 
- . NAAN. OE 
WES UN NOIINONININN GO 
tetyvisrestertsr 
Ce 


GgourCEe CC oOCCOCOOC Cor 


GO 
ae) 
ae 

Ta) 
Ce 


CeG2CSURD OC 
“0 .4809C8D GO 
0 

co 

a0 

-C.4809C8D OC 


0 64809C8D 
= .4ARU9L PD 


oo 
fn 


SRO9UHRD 

90 282D 

TPO 
420999D OC 

4BI9D9RD 

+P 0908) 

9.489998) 
ce 


—C 248i 


a) 


c 
—§ 4 8C9UED 
-C.6AC9NRD 20 


~U eG ANG 
-0 4809099 
0 648099 20 


-0 24809080 


=C 


Et tt tt tt tn 
CEVCOCCOCc BCOOCOOT COC 


eqsonesasaco0eosacea: 
SR RRRRERERR SR ERR 
atreswrerysersrrtssr se. 
tereereterrsrtsese 
cakeniemien hen ienkanhethanhan ten hanhstien len lan kanla 
lopelolalolaleelolale! alelelalajlel 
Et tt et et ad Net tet tr 
eter er eee eee eeeeens 
cocoocoso0cococe oc o« 


OoorNe LHMNNOOSS 
mt tt et tt tt cmt mel en meh et 
CCR TE PERT EST Ee 
oOs35306uq0 Ss30900650ca 
AEF ALZONF MOWCON SD 
FAKROTH CANVOMeEmOM 
DOROOD LHOMO HAWS 
Re vom woReNntaw cae 
CODORD® etry tu comm 
OF SDH INO RO NNR Ome 
eeteeveeveeetr eee eeteeee 
ce a 


soceescsecs: 
baa i ira 0h site 


et a td tt pk ot ah a 


eocooocoooccuo0ococt 


scooococoacaasacaaac 
NCNNOMNEOMD DOO OHOwe 
ONG MHOREREREESREREARS 
Volos atl alaralal a al al al al al al al an 
DIF OODDNDODOHOMOMMOMaMc 
OR RRRREERARRARERERAES 
PITT SSS STS sstss. 
eeeeeeeeeeeeee eee 
lelelelelelelelelelelelelelelelelald 


ANM SHOE DROANM SHO C 
te ime 


Se 


te tt 


owe Ge Sree eo 


ae 
=O G59998D 


> .4299780 


=f .4anenan 
=f .42999cRD 


A tt tt 


"Cc COoucceouw 


396005959000 
= REAR RARE 
rrt@rrrerse 
res—r~rwrs rs 
tt et mt te 
IOC oedIcacs 
tt tt et mt et 
. a CI ee 


, 
) 


-00740149U-1 
-022590520- 


tt et at tt 


del alejolqlolai« 


IGBO060C00 
LoS MUBeOrOD 
Seo eel eel Ded eel ed el el 
lal sales al al ata 
OMnmmarnen 
le 
taerrrrts 
een eee eee 


scaoo0o0co 


SCOANMT HO 
A RINNINOINNN 


COEFFICTENTS 


FOURTER 


AXSRES 


, 
eeoereree te he eoeeseoes 
IOS 2 IVECO 


oeeeer et oe 


O@usco Cue - wrt ou 


God. € ce yure > e¢ Cee SOS 
eee ee eee eee ee eee were eres ees 


OOOSNCCOaAIVGC TOVCOCE VOGSCCOCCS 


Cw ae © SUSE SOY Och t wot oer 
SSEE SS Se CaS eee Sere See eS 


COOL OTe WC OC COLCOrEocnoece 


at OU PHO WU rs UHL NHOOY ANY 
Et et et me et et et 

DPECECERE CERES SRE Cee EE GE 
2 CH0S5D0C00a6cag5S06G06 


mTOMDenmrNOIWmMOLS 
OWWEWDENOOKRDPOOO 


0.881757D 
te 7638330 


OOO9OCS” OVUCCCOOOCSCEr COOTOOG 
eeereteereeereeeeeeeeeeeeeeeeve 
MOCSCOOGCIOCE COCOYVLEGEOSOSCCOCH 


SOOCCOCOOCT COOOOCTC COOOCCCOGSO 
eet eet eee eee eee eee essen eee 


CoocoocoooooooCoCoCcocoo00ecocono 


ANM SHOR DRORNM SHORE OOMNE THO 
tt me tt ONIN 
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ARADSY 


5 ad acl ad et Oe Ol Dd 


et 


eaae level 
ert (8 G1L0=-10 

20-16 
ec OIl¥4sUu-hl 


“eSTLeeou-lt 


eB AR Wey 


43452 50-17 


THT LooOuUK17 
& 


le 


} 


eva (299017 
ech PV suHt 


JIesNTITIIsNv ITI sosesrss 
Lo FARA AA HR Rt A te tte 


malay 


=) AGOAOCSCGIS 
am fy oA Ne mas LrC 1: 
Chwoeaer mre 


\OONC Mure 
OTAUMBHF Her 
aM mOaproce 
OMG OND ~ 


6640 


pS 
ce 

hu 
«7 
cs 


621914020 
624911030 
te aay 


—tie2t 


b-07 
IG 3RaD—-C8 


2204 


iy oa iar Ari 
PGK ' 


fara) 


4p 
OeLAM% ZOU HA14 


Dae DOR MN IMM O- 
CaO GIGS OUNCE Cd ee aes 
NA TONKR ta OT Tm DIN 
PONTO OR MND AH OMMNOE OS DC 
DOOM AS DR RNA ENDER See 
eeeeeveeenereeteeeeeeee ee 
BOCOCOCCOGC SICOSSUOGG 
at 1 | a pit 


04809; 


af 


AOD KANN CHOMY HDWMWNNMW on NW 
Pa eae | 
eogqgcesq0csg0cceog005g0G¢ 
NANO ROW OTNHORDATL GR AORDS 
DOME MAOMOTS DMOANKOMANH- 
BMPACC LSS OLOG PA QHOKNHOMM. 
Baw Ou ODM CHUNG SENN 
WO MMRRROR DOR YeNwOI NOt 
MOR AOL ROO NM Dm OM HORNE 
oe ee ee ee oe ee ee eee 
set ah ah CAS oe 


eh eae hs 


5 
2 


eral al aralal at oral al are) Sioa al al am ol al als 
EE el pk tt tt 


6239920-1 
N.327840U-16 
pO ten ha 


“02199247015 


-0 


NNR NAN 
Titiiia. 
i i 


RNNRUNNONNN 


60- 


OtTaeKrno-c 

NORor oon 
NN Se eter 
eeeeete 


TrPeTT 


MOKROKNMSNOFDOOMNM Su 
tt et INNIS 


-0.2382720-02 


-0. 


CG 


bak © 


=" 


i 
c 


Yes 
Ps 
rs 


— 4+ 


2141-15 


+ IO 


weak, ED? 


We ae BOE ve 
“0. 79152R5-95 
WALL DATA aT 


’ 


6 


“TATION 


Ny tt tt mt et ter 


Xe 355271 0-1¢ 
Oe I562TlU-le 


3uR RE af 
ae 3 yteayy 
F436 Aly 
aurrmsay 
44,476 ary 


9D 


COCOCoSe eres 
SOGCSC ee OSSLP-~ 
COCOCYOOCeE OO 


19UCOCD 


LONC 
1a COCOOH 


EN th ht tt tt tN det od tO Oe tt 
LO DOGVPOCOVSESCOCELAGOCOTCGOIOCS 
EEG ERELL PECS NED EP Ee Tee ee 
SGSOEADCEHOSCADOOOCOAGHaGCoOCHAa 
SE ed ed et nt dE gee ed ed ed eh et ed emt ed md ed eet tl et 
FFIFRMF OV FROCICCSHSCCAOSISSCOTZO 
PEO OSL OOS CSS Ooo S S 
ewe SOO wee a 


"0 
OO 
cc 
ul 
uc 


theta 


REDD. Ga 
CO2CCS 


30450 
3045) 0 


61404650 


DHOHOOH GO QLL OO OLLH OOOH OOL OOH 
ee 
CUCSCOeS DOCCOCUCOC OOCUULOC OC 


£661 27045L 


O.6170650 
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